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Abstract 
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1 Introduction 



After dilute condensates of bosonic atoms were produced and observed in the laboratory 
by Cornell considerable interest arose concerning the dynamics of the more complex 
system formed by two coexisting, coupled condensates H, El- I n Te ^- 11 the dynamics 
of component separation of a magnetically trapped dual condensate has been studied by 
making use of the possibility of adjusting independently the trapping conditions for each of 
the two components. In this case the composition of the mixture remains fixed, and the 
finite size of the system plays an essential role. Other situations, in which the coupling 
allows for interconversion between the two condensate types, have also been considered from 
a theoretical point of view. Here, even under equilibrium conditions for which the spatial 
dynamics of the trapped condensate particles is essentially frozen, one still has to consider 
the interesting collective dynamics of the composition of the dual condensate. One of the 
situations that has been considered in this connection involves two different internal states 
of the atoms, as in refs. |2|, [|. In this case, the interconversion coupling is provided by laser- 
induced Raman transitions between these two states 0. More recently, the observation 
of effects of the so called Feshbach resonances in atomic Bose condensates H led to the 
consideration of possible experimental situations involving coupled atomic and molecular 
condensates, the latter occurring in the two-atom channel responsible for the resonance 
phenomenon || [7J. In this paper we explore the fact that, by making use of the assumption 
of frozen spatial dynamics, the dynamics of condensate composition can in both cases be 
treated in terms of simple collective variables which evolve in an essentially classical regime, 
quantum numbers being of the order of the number of atoms involved in the condensate. We 
can thus derive classical, canonical equations of motion governing the model composition 
dynamics. 



2 Effective dynamics of coupled condensates 

Although from a microscopic point of view the dual trapped condensate dynamics involves 
many quite subtle questions of atomic physics, once these are duly tamed to the point where 
they can be manipulated in the laboratory it is possible to encapsulate their effect in a few 
dynamical parameters for the purpose of studying the overall behavior of the condensate. 
Theoretically, the two above mentioned cases have thus been modeled in terms of an effective 
Hamiltonian density of the form 0, |6], 0] 



7i = 7~l a + l~Lb + T~C C 
where the first two terms describe the individual condensates, i.e 

^ 2 V 2 X a ,b , , 
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and the last term contains the coupling of the two condensates. For the atom-atom (AA) 
case it is written as 



h[ aa) = \ \<f) a n <pb r +a(<p*j b + <p*m (3) 

while for the atom-molecule (Feshbach resonance, FR) case it reads 

nf R) = \ \ 4> a ri h r +«(c 2 06 + <pi<pd- (4) 

The parameters e a ^ represent possibly different intrinsic energies in the two boson channels. 
Their elastic interaction is described in terms of pseudopotential parameters \ a ,b, A and 
channel coupling is represented by the parameter a. In the AA case considered by Zoller || 
the latter involves the intensity of the laser responsible for the Raman transitions, and should 
therefore be considered as an externally determined control parameter. In the FR case, on 
the other hand, a stands for the coupling to the quasi-bound molecular state responsible for 
the resonance and therefore determining its width, while e& — 2e a represents the "detuning" 
away from the resonance, controlled externally by means of an applied magnetic field || . 

For sufficiently extended and uniform systems the effective field operators (f> a , <pb can be 
usefully expanded in a momentum basis as 



k v k 

where the nonhermitean mode operators oj, satisfy standard bosonic commutation rela- 
tions. Furthermore, the depletion due to correlations being small in low density systems ||, 
in the condensate regime essentially all the bosons are in the zero momentum mode, so that 
the relevant Hamiltonian density reduces to the simple single-mode form (for simplicity we 
omit the zero momentum label in this case) 

H a — ► y a a + T^^a'a'aa, 

a similar expression for Tib and 

Hi AA) -> ^a^ba + ^(a^b + tfa), (5) 

H{ FR) -> Aatfttfta + " ( a t a t 6 + 6 t aa ) ( 6 ) 

for the atom-atom and Feshbach resonance cases respectively. In both cases the total number 
of atoms (namely, a^a + b^b in the AA case and a) a + 2Wb in the FR case) is obviously con- 
served. Exact single mode coupled condensate solutions can therefore be obtained through 
the diagonalization of finite matrices. Furthermore, the single mode coupled condensate 
Hamiltonians belong to the special class of the so called Curie- Weiss models M, for which 
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the mean-field approximation becomes exact in the thermodynamic limit (number of parti- 
cles — > oo at constant density). In order to explore this latter feature, it is convenient to 
express 7i^ AA ^ and 7i^ FR ^ in terms of alternate dynamic variables, which can be conveniently 
chosen so as to obey standard SU(2) commutation relations. 



2.1 Two atomic condensates 

In the AA case, the appropriate variables are just the well known Schwinger realization of 
the SU(2) algebra in terms of two types of bosons [[UJ 

J z = ^(a f a - b ] b)- J + = A = a t&. 

In this case the Casimir operator J 2 appears as J(J + 1) with J = h(a'a + b'b), half the 
number of atoms, and the Hamiltonian density can be written as 



= e ^j+ e -^j z + ^- 2 (j + j 2 )(j + j z -i) + ^- 2 (j-j z )(j-j z -i) 

In order to take the thermodynamic limit it is more convenient to work with the scaled 
Hamiltonian = VHS AA) / J which reads 



h (AA) 



e a + e fe + (e a - e b )— + ( 



J, 



J, 



(7) 




where the atom density n = 2J/V has been introduced. The thermodynamic limit consists 
now in letting J — > oo at constant n. The spectrum of the scaled components J-J J remains 
bounded in the closed interval [-1,+1] and becomes increasingly dense as J is increased. The 
thermodynamic limit corresponds therefore to the classical limit of a dimensionless angular 
momentum-like algebra which may be formally characterized as 



lim — 

J^oo % 



Jk 
1' 



— \jkijl\ ^klmjm — 



lim 



J 71 



J-+oo J 



where the curly brackets now denote a Poisson bracket. In the thermodynamic limit, Eq. 
(0) yields the (quasi-) classical expression for the scaled energy 
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nXn . nA 



4 v • 4 



T? A 

+ Y (l-f z )+2a 3x . (9) 

Alternatively, one can express Eq. (|9]) as a proper Hamiltonian function, in terms of a pair of 
canonically conjugate variables. One such pair which is quite convenient has been previously 
obtained using the quantum kinematical scheme developed by Schwinger in terms of unitary 
operator bases |TT| and the associate (discrete) Weyl-Wigner transforms [jnj . After carrying 
out the thermodynamic limit in the way just described, one finds that the (angle) variable 
which is canonically conjugate to the ("action") variable j z is just the azimuthal angle tp, so 
that j x can be expressed canonically as 



jx = V 1 _ Jz cosl P- ( 10 ) 

Using this fact, one may derive from Eq. (|]) equations of motion for the variables j z and 
(f. A technical point to be observed here is that the dimensionless character of the "action" 
variable j z corresponds to measuring time in inverse energy units. In order to introduce an 
appropriate time scale one may write, consistently with Eq. (M), 



dj z .. d J z J 
— — = hm — = hm — 

dt j^oo at J j->oo ih 



j ' j 



h {AA) 

a>^-> en) 



which gives 



dip dh^ A) e a -e b n\ a n\ b n\ . 2aj z 

+ -^-(1+ Jz) - -^-(1- Jz) - ~rJz / cos y 



dt hdj z h 2h ' 2h ' h j- 

(12) 

dj z dh%*> _2a h _ 2 



, x — K sin tp. 

dt hdip h y Jz Y 

In order to explore further the significance of the canonical variable ip it is useful to derive 
the second Eq. ([12]) in the following alternate way. The (zero momentum) field operators 
4> a and (f) b satisfy the coupled nonlinear equations 



iU(f) a = (e a + X a | <p a | 2 +A | 4> b | 2 ) <p a + a<j) b 

iftipb = (e b + X b | (p b | 2 +A | (p a | 2 ) (p b + OL(p a - 



From these one may easily derive the equations for | <j) a and | <p b 
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i%^2 | 0a | 2 = 2za Im0 6 0* 

(13) 

^rf£ I 06 I 2 = ~ 2ia lm 0b0a 

which show an explicit dependence on the phase of 0b0*. Writing this object as | 06 1 1 a | 
e tS(AA) the difference of these equations appears as 

^(|0 a | 2 -|0j 2 )=4a|0j|0 a |sin^) 

while their sum merely gives the conservation of the total number of atoms. Using the 
definitions of J and J? one has 



I r u I I tv | — j i | t u || f u, | — 

so that the difference equation becomes 



2 ~ 2Jz - 1 ' 1 ' ' / (J + ^)(^-^) 



This corresponds to the second Eq. (|T2"D and furthermore identifies the angle variable with 
the phase 5 ( -" 4A - ) . 

The density dependence of the Hamiltonian Eq. (^) also gives a direct analytical expres- 
sion for the pressure exerted by the interacting condensates. One has in fact 



p(AA) _ n 2 



dhi^ n 2 



qc 



A a (l + j,) 2 + A (l - ]z f + 2A(1 - j 2 z ) (14) 



9n 4 

In order to avoid the collapse of stationary states this expression must be positive when 
evaluated at the equilibrium value of j z . 



2.2 Feshbach-resonant atomic- molecular condensates 

The FR case, on the other hand, can be treated in precisely the same way once an appropriate 
realization of the SU(2) algebra is constructed. To this effect, consider the case in which 
there are N atoms present. The relevant finite dimensional space in which ri^ FR ^ is to be 
diagonalized is then generated by base vectors of the form | N—2nb, rib), where the two labels 
denote the number of atoms and the number of molecules respectively. For definiteness, iV 
will be assumed to be even, so that runs from zero to N/2. The next step is to identify this 
basis with the J z eigenstates of a SU(2) multiplet associated with the eigenvalue N/4(N/4+l) 
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of the Casimir operator. A convenient way of doing so is to identify the eigenvalue of J z 
with n b — N/4, so that J z = \(2b% — a^a) and 

N N 
| N - 2n b , n b ) <->| J = — , J z =n b - —). 

If one then defines J± in terms of their standard action on the J z eigenstates, i.e. 



J± I J, Jz) = \JJ(J + 1) - Jz(Jz ±1)\J,J Z ± 1) 

one finds, after a straightforward calculation, 

1 1 

J+ = ^=ggtf; J_ = a'a'b- 



'2(tfa + 1) v /2(ata + 1) 

The FR Hamiltonian density can now be written in the form 

-H (m = ^^J + ^^J, + ^2(J-J,)(J-J 2 -1) + ^(J + J Z )(J + J Z -1) 

+ ^2(J 2 - J 2 ) + ^ J - J z ) + 2 + \J4{J -J z )+2 J, 

and the corresponding J-scaled Hamiltonian, after taking the thermodynamic limit, becomes 



h%*> = lim^V- (15) 



= e b + 2e a + (e 6 - 2e a )j z + -^(1 - j,) 2 + -^(1 + j 2 ) 2 + — (1 - j 
+2aty/ny/l _ ^ ^ 

where the vanishing of the commutators [^f, i n this limit has been used to obtain the 
last term. An interesting feature of these Hamiltonians is the quenching of the coupling 
term involving a when the molecular component is dominant (J z close to J or j z close to 
1). It can be traced ultimately to the feature of Bose statistics which associates the factor 



(N — 2n b + 1)(N — 2n b + 2)n b to the amplitude for converting one molecule into two atoms 
when there are n b molecules present. This factor in fact decreases faster than that associated 
with the plain J_ operator when n b is large (close to its maximum value N/2). 

One may next adopt again the canonical form of j x given in Eq. ([10]) and the appropriate 
time scale to write the equations of motion 

dip e b -2e a n\ a . nX b . nX . ay/n 1 + 3j z 

Tt = — n {l ~ 3z) + -4h {1+3z) -T 3z -—VTTTz COSlf 

(16) 

dj z 2ay/n 

dt h 1 



1 -jz)\l + Jz sin if. 
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In a way which is completely analogous to the AA case one may here identify the canonical 
variable ip with the phase 8^ FRS > defined as 



U 2 |«M e**«. 



The pressure can also be obtained in this case with the result 



p(FR) 



n 

T 



Aa(l-J 2 ) 2 + ^(1+J 2 ) 2 + A(l-J 2 2 



+ an 



3/2 fl _ 



J.: 



J z COS 



(17) 



which reproduces (with different notation) the result obtained in ref. [[]]. When the square 
brackets are positive for the equilibrium value of j z , one may still have a domain of negative 
pressures at low densities when the coefficient of n 3//2 is negative. Also as observed in ||, 
this implies that the system becomes "self-bound" and saturates at the value of n for which 
the pressure vanishes. 

In order to make the model sufficiently realistic for the FR system it is important to take 
into account the loss of atoms which may be expected due to an enhanced rate of inelastic 
collisions involving the molecular channel directly |5j, [7J. This loss has in fact been used as 
a signal to detect the Feshbach resonance in ref. ||. The simplest way of including the loss 
of atoms consists in adding appropriate master loss-terms to the equations of motion for 
the atomic and molecular densities, n a and n\,. Under loss-free conditions, we may use the 
expressions 



n„ 



1-3* 



rib 



-n 



and the equation of motion for j z , Eq. (|T^), to write 



dn a 
dt 
dnb 
~~dT 



ndj z 
2 dt 

ndj z _ 

4 dt 



an 



3/2 



h 



(! - 3z)\fl+ jz simp 



an 



3/2 



2h 



1 - j z ) V 1 + Jz sinv?. 



These equations are the FR analogs of Eqs. (|T3D, used to identify the canonical variable 
ip with the phase of an appropriate combination of field variables. To take loss effects into 
account these equations are replaced by 



dn a 
dt 
dn h 
~~dT 



an 



3/2 



an 



h 

3/2 



1 - jz) V 1 + jz sin if - c aa n a - c ab n b n a 



2h 



1 - j z ) Jl + j z sin (p - c bb n b - c ba n a n b 



8 



which involve four loss-rate coefficients Cjj. These are defined so that CijTij represents the 
decay constant of the z-boson density due to collisions with j-bosons. The important decay 
rates in the molecular channel thus imply larger values of and Q, a when compared with 
c aa and c a b. The time evolution of the total density n is now determined by the equation 

ii = 2 ir + ir = - ! j I 2Caa(1 - 3zf + Cbb{l + 3zf + {Cab + Cba){l ~ 3l) \ (19) 

which must be solved together with Eqs. flUf ) in order to include atomic losses in the 
dynamics of the interacting condensates. It should be noted that this procedure corresponds 
exactly to that adopted in ref. M, where one works with the nonlinear equations of motion 
for the fields (j) a , <pb modified to have complex pseudopotential parameters. In fact, the terms 
involving the imaginary parts of these parameters cancel from the equation of motion for 
the variable j z , but give non vanishing contributions to dn/dt which exactly reproduce Eq. 

ED- 



3 Numerical results 

In order to explore the classical propensities of the quantum, Curie- Weiss coupled condensate 
dynamics, it is most convenient to make use of a phase-space quantum description. An 
appropriate description of this sort, for systems evolving in quantum phase spaces of finite 
dimensionality, consists of the discrete action-angle Weyl-Wigner representation used in Ref. 
fT2|| . In the next few lines we simply collect the final prescription to obtain the discrete 
transforms, and refer the interested reader to this reference for further details. For a state 
| a) represented in the appropriate multiplet base | J, m) as 



a) 



j 

E 

m=—J 



a„ 



J, m) 



one first constructs the matrix 



V ZJ + 1 m=-J 



2ni I N 

k(m H — , 

2J+1 v 2' 



where the range of the integers k and / is — J < k,l < J and the index {m + /} denotes the 
value of m + / cyclically confined to the range — J, J of the basis labels. Explicitly, one has 



{m + 1} = m + I - (2 J + 1) Floor 



'm + I + J" 
2J+1 . 



where Floor(x) denotes the larger integer (negative for x < 0) less than or equal to x. The 
desired discrete Wigner phase-space representative a w {p, q) of the state | a) is then obtained 
as the double (discrete) Fourier transform 
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(2 J + 1) a w (p, q) = ex P 
v 2J + 1 k,l 

In this expression the range of the integers p and q is also bounded as — J < p, q < J, 
and the properly scaled variables corresponding to j z and if are q/J and 2np/{2J + 1) 
respectively. The Weyl transform of the Hamiltonian can be obtained in exactly the same 
way, replacing the amplitude products a m a| m+/ j. by the matrix elements (J,m \ VH./J \ 
J, {m + I}) multiplied by the number of states 2 J + 1, when evaluating r(k, I). 



2m 



2J 



-(pk + ql) 



r(k,l). 



3.1 Two atomic condensates 

We restrict our numerical treatment of the composition dynamics of two atomic condensates 
to the "symmetric" case, in the sense of ref. i.e. e a = = (implying a suitable 
definition of the energy scale) and A a = Xb = Xq. Furthermore, we adopt units such that the 
atom density n = 1 and H = Xq = 1 so that energies are given in units of Aon, times are given 
in units of H/XqU and the various possible situations will unfold by varying the remaining 
parameters, A, a and, when short of the thermodynamic limit, J = N/2, cf. Eqs. (0) and 

(!)• 

The numerical spectrum of h^ AA \ Eq. ([?]), is given for A = 1.5, a = .1 and J = 
20 in Fig. 1, together with the mean b- type atom numbers (n^) = (tfb) and variances 
a b = J ((tfb) 2 ) — (n b ) 2 for each of the corresponding eigenstates. The states with energy 
Ek < 1 are in fact nearly degenerate doublets (see Fig. 5 below) of "Schrodinger cat states" , 
consisting of superpositions of basis states with (rib) ~ and (rib) ~ N 0]. The large values 
of erf, for these states signal the large, strongly correlated fluctuations of rib an d n a . The 
doublet structure can be immediately understood with reference to the Weyl transform of 
hS AA \ shown in Fig. 2 together with the quasi-classical energy surface corresponding to h^ AA) : 
it results from the symmetric minima near j z = ±1 separated by the barrier which peaks 
at j z = 0. As also discussed in ref. M, this situation results from having A > Ao and small 
enough a ("weak laser" case). The energy splitting of the doublet members approaches zero 
as J is increased, leading to degeneracy in the thermodynamic limit. The Wigner functions 
corresponding to each state in the lowest doublet for J = 20 are shown in Fig. 3. The 
lowest (highest) member of the doublet involves a symmetric (anti-symmetric) superposition 
of states | J, M) strongly concentrated on the largest values of | M | . The Wigner functions 
show moreover that these states are strongly peaked also in the conjugate angle variable, as 
a result of the dips of the energy surface at if = ±7r. Increasing a will depress the pass along 
the (p — 7r line leading eventually to a minimum at j z = and therefore to ground states 
dominated by small | M | components. This is shown in Fig. 4, where the probabilities 
p(M) =| ( J, M | g.s.) | 2 are plotted against M for a range of values of a. 

A most dramatic consequence of the strong collective character of the composition dy- 
namics as described by the Hamiltonian h^ AA \ Eq. ([7]), is that, as was pointed out before 
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(see e.g. ref. [[Uj), a suitably scaled distribution of the energy eigenvalues quickly approaches 
the quasi-classical distribution in which the fraction of levels with eigenvalue smaller than 
a given energy is proportional to the area of the phase-space domain in which the quasi- 
classical energy is less than this value. This ultimately allows associating energy eigenstates 
with phase-space trajectories on an essentially one to one basis. Fig. 5 illustrates this feature 
showing the fraction of the total number of states having energy less than E as a function 
of E for several values of J together with the quasi-classical limit. One sees there that even 
the case N=20 is already fairly close to the quasi-classical limit from which it differs mainly 
by an overall translation in energy related to the terms of order 1/ J in Eq. (|7|) . Also visible 
is the staggering associated with the doublet structure of the A > Ao, "weak laser" case, 
which persists up to energies close to that corresponding to the separatrix going through the 
highest point of the ip = tc pass at j z — (E — 1.05 for this example, in the quasi-classical 
limit). An enhancement of the level density at this energy is also clearly visible. 

3.2 Feshbach-resonant atomic-molecular condensates 

While in the AA case the relevant external control parameter is the intercondensate coupling 
a, associated with laser intensity, in the case of realistic FR systems the relevant parameter 
is the detuning 5 = e a — 26b, controlled by an externally applied uniform magnetic field. 
In order to illustrate the composition dynamics in this case we take A a = A^ = A = 1 and 
again choose units such that h — n — 1, the zero of the energy scale being defined so that 
e a = 0. The remaining parameters to be considered are therefore a, here related to the 
resonance width, the detuning parameter e?,, which is the intrinsic energy difference between 
one molecule and two atoms, and (again when short of the thermodynamic limit) J = N/4. 

The numerical quantum spectrum of VHS FR > / J is given in Fig. 6 (similar to Fig. 1, 
corresponding to the AA case) for a = \[2/2, = 2 and J = 20. The mean molecule 
number (rib) and variance o\> for the eigenstates are shown at the corresponding energy 
eigenvalues. The Weyl transform of the quantum Hamiltonian is shown in Fig. 7, together 
with its classical limit, Eq. (|1^). The scaling property of the level density with J, analogous 
to that shown in Fig. 5 for the AA case, is shown in Fig. 8. 

In this case, the quantum ground state is strongly localized, both in action and angle 
variables, in the non-degenerate minimum of the energy surface, as shown by its Wigner 
function, Fig. 9. Increasing the value of the detuning parameter energetically favors the 
lower M components in the ground state, which becomes eventually an essentially pure 
atomic condensate, as can be seen in Fig. 10, where = 50 with no change in the remaining 
parameters. This figure shows the Weyl energy surface for this case and also the Wigner 
function of the corresponding ground state. Since this energy surface is dominated by the 
detuning term, and therefore rendered only very weakly dependent on the angle variable, its 
ground state appears as strongly delocalized in this variable for J = 20. The definition of 
this ground state in the angle variable is however expected to increase for larger values of 
J, as a more collective state fits in the very shallow minimum which persists at ip = n, near 
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jz = "I. 

Next we use the e& = 50 quantum ground state as an initial condition in order to obtain 
numerically the (loss free) quantum evolution under the e b = 2 Hamiltonian. This corre- 
sponds to the scheme devised in ref. J/J, in which the ground state is prepared at some 
initial large value of the detuning which, is then suddenly reduced by adjusting the applied 
external magnetic field. The absolute square coefficients of the decomposition of the initial 
state in the e = 2 eigenstates are shown in Fig. 11, and the resulting time dependence of 
the mean value of J z j J for J = 20 are shown in Fig. 12 together with the corresponding 
time-dependent variances. The mean number of atoms and of molecules, (n a ) and (rib), are 
related to (J z )/J simply as 



(n a ) = 2J 



1 - 



J 



(n b ) = J 



1 + 



J 



and will therefore oscillate with opposite phases as a result of the conservation of the total 
number of atoms in this calculation. Also shown in the same plot is the function j z (t) 
obtained by integrating numerically Eqs. (pT6|) with e b = 2, a = y/2/2 and using the phase- 
space coordinates of the quasi-classical equilibrium point corresponding to 5 = 50 as initial 
conditions. Note that, in view of Eq. (|TT|), the appropriate time parameter to be used in 
the quantum calculation in order to compare with the quasi-classical result corresponds to 
solving the Heisenberg equation of motion 

J z VHSFR)' 



ih d 


Jz 




~Jdi 


. J _ 





J 



with J = 20 in this case. The similarity of the two graphs illustrates the relevance of the 
quasi-classical calculation of collective properties of the system already for J = 20. 

Finally, we include the phenomenological loss terms in the quasi-classical equations of 
motion by integrating simultaneously Eqs. (|T6|) and Eq. fll9|) with the realistic values of 
the loss-rate coefficients used in ref. @, namely c aa = c ab = and c bb = c ba = 0.5, in the 
units adopted here. Fig. 13 shows the resulting time evolution of j z {t) and of n(t), the 
latter quantity being measured in units of the initial density n(t = 0) = 1. The atomic and 
molecular densities are now given by Eq. (|T8|). Their decrease is due mainly to the decrease 
of the total density n since, as can be seen in Fig. 13, the amplitude of the oscillations of j z 
remains relatively immune to the loss effects. Furthermore, the undulating behavior in the 
decrease of n(t) reflects the asymmetry of the adopted values of the loss-rate coefficients, as 
a result of which the losses take place in the molecular channel only. This in fact causes the 
effective loss rate of n(t) to vanish as j z (t) approaches its lower bound. 



4 Concluding comments 

The fact that the composition dynamics of essentially undepleted, coupled uniform Bose con- 
densates can be completely reduced to a very compact form in terms of collective variables 
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only is a direct consequence of the Curie- Weiss character of the underlying model Hamilto- 
nians. While the terms included in these Hamiltonians have been restricted here to those 
considered earlier ||, the inclusion of other terms, such as two molecule collisions leading 
to four condensate atoms and atom-molecule collisions leading to three condensate atoms in 
the Feshbach resonance case, is completely straightforward as they can be readily expressed 
also in terms of the same collective variables. Not only the resulting collective Hamiltonians 
reduce to classical form in the thermodynamic limit, but it it also an easy matter to obtain 
complete numerical quantum solutions short enough of this limit. The study of these solu- 
tions indicates that the quasi-classical thermodynamic limit is approached rather fast, while 
allowing for insight into the dynamically relevant correlation processes. 

It should be kept in mind, however, that purely collective dynamical treatments of this 
kind involve a severe truncation of the full model problem as stated e.g. in Eqs. (|l|) to (f|). In 
particular, collective excitations will in general couple to depletion effects which fall beyond 
a purely collective description. We believe nevertheless that the remarkable simplicity of the 
latter readily provides useful information on such systems and may give useful indications 
for the development of more complete treatments. 
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Figure Captions 



Figure 1. Energy spectrum of hS AA \ Eq. (0), for a symmetric, "weak-laser" case. Units and 
parameter values are explained in the text. The points below E ~ 1 correspond to nearly 
degenerate doublets (cf. Fig. 5). Mean numbers of b atoms and the respective variances are 
shown at the corresponding energy eigenvalues. 

Figure 2. Discrete action- angle Weyl transform of the Hamiltonian used in Fig. 1, (a), and 
the corresponding quasi-classical energy surface, (b). 

Figure 3. Discrete action-angle Wigner functions corresponding to the lower (a) and upper 
(b) member of the lowest doublet in Fig. 1. 

Figure 4. Evolution of the probabilities p(M) of the various | J, M) components of the 
ground-state as a function of the "laser strength" control parameter a. The double-peaked 
distributions correspond to "Schrodinger cat states" . 

Figure 5. Evolution of the integrated level density as a function of J. Other parameters 
are as in Fig. 1. The quasi-classical limit is also shown. 

Figure 6. Typical "small detuning" energy spectrum of V7i^ FR ^ / J. Units and parameter 
values are explained in the text. Mean numbers of molecules and the respective variances 
are shown at the corresponding energy eigenvalues. 

Figure 7. Discrete action-angle Weyl transform of the Hamiltonian used in Fig. 6, (a), and 
the corresponding quasi-classical energy surface, (b). 

Figure 8. Evolution of the integrated level density as a function of J. Other parameters 
are as in Fig. 6. 

Figure 9. Discrete action- angle Wigner function corresponding to the ground state of the 
spectrum shown in Fig. 6. 

Figure 10. Weyl transform of the Hamiltonian, (a), and Wigner function of the ground 
state, (b), for "large detuning" e& = 50. Other parameters as in Fig. 6. 

Figure 11. Squared expansion amplitudes of the "large detuning" ground state of Fig. 10 
in terms of the "small detuning" eigenstates corresponding to Fig. 6. Results are given for 
the full J = 20 spectrum. 

Figure 12. Exact quantum evolution of the "large detuning" ground state under the "small 
detuning" Hamiltonian for J — 20. The heavy curve shows the time evolution of (j z ), and the 
error band shows the time evolution of the root-mean-square dispersion of j z . Also shown 
(white squares) is the result of an integration of the quasi-classical Eqs. flIB|). See text for 
details. 

Figure 13. Numerical solution of the coupled quasi-classical equations flTE| ) and (^). See 
text for details. 
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